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Abstract 



Using the trace anomaly relation, low-energy theorem and Witten-Veneziano 
formula, we have developed an analytical formalism which allows one to cal- 
culate the gluon condensate, the topological susceptibility and the mass of 
the 7/ meson in the chiral limit as functions of the non-perturbative vacuum 
energy density. It is used for numerical evaluation of the chiral QCD topology 
within the QCD vacuum model consisting mainly of the quantum component 
given by the recently proposed zero modes enhancement (ZME) model and 
the classical component given by the random instanton liquid model (RILM) . 
We sum up both contributions into the total, non-perturbative vacuum en- 
ergy density. A very good agreement with the phenomenological values of the 
topological susceptibility, the mass of the 7/ meson in the chiral limit and the 
gluon condensate has been obtained. This puts the above mentioned QCD 
vacuum model on a firm phenomenological ground. 

PACS numbers: 11.15.Tk, 12.38.Lg 
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I. INTRODUCTION 



The nonperturbative QCD vacuum has a very rich dynamical and topological structure 
[1]. It is a very complicated medium and its dynamical and topological complexity means 
that its structure can be organized at various levels (quantum, classical) and it can contain 
many different components and ingredients which may contribute to the vacuum energy 
density, the one of main characteristics of the QCD ground state. The quantum part of 
the vacuum energy density is determined by the effective potential approach for composite 
operators introduced by Cornwall, Jackiw and Tomboulis (CJT) [2] (see also Ref. [3]). It 
allows us to investigate the non-perturbative QCD vacuum, since in the absence of external 
sources the effective potential is nothing but the vacuum energy density. It gives the vacuum 
energy density in the form of the loop expansion where the number of the vacuum loops 
(consisting of the confining quarks and non-perturbative gluons properly regularized with 
the help of ghosts) is equal to the power of the Plank constant, h. 

In our previous work [4] we have formulated a new, quantum model of the QCD ground 
state (its non-perturbative vacuum), the so-called zero modes enhancement (ZME) model. 
It is based on the existence and importance of such kind of the non-perturbative, topo- 
logically nontrivial quantum excitations of the gluon field configurations (due to the self- 
interactions of massless gluons only, i. e. without any extra degrees of freedom) which can 
be effectively, correctly described by the g _4 -type behaviour of the full gluon propagator 
in the deep infrared domain. The correct treatment of such strong singularity by the dy- 
mensional regularization method within the distribution theory was one of the highlights 
of the above mentioned preliminary publication [4]. Our model allows one to calculate the 
non-perturbative vacuum energy density from first principles using the CJT approach for 
composite operators [2]. We have also formulated the method of how to determine numer- 
ically the finite part of the vacuum energy density. We propose to minimize the effective 
potential at a fixed scale as function of a parameter which has a clear physical meaning. 
When it is zero then only the perturbative phase remains in our model. Equivalently one 
can minimize the corresponding auxiliary effective potential as a function of the ultraviolet 
(UV) cut-off itself. The non-perturbative chiral QCD vacuum is found stable since its main 
characteristic -the vacuum energy density - has no imaginary part and it is always negative. 

Within the ZME quantum model of the QCD ground state [4], the vacuum energy 
density depends on a scale at which the non-perturbative effects become important. If 
QCD itself is confining theory, such a characteristic scale should certainly exist. The quark 
part of the vacuum energy density depends in addition on the constant of integration of the 
corresponding Schwinger-Dyson (SD) equation. The numerical value of the non-perturbative 
scale as well as the above mentioned constant of integration is obtained from the bounds 



for the pion decay constant in the chiral limit by implementing a physically well-motivated 
scale-setting scheme [4]. We have obtained the following numerical results for the non- 
perturbative vacuum energy density, e = e g + Nf€ g , 



87.2 < F° < 93.3 (Mel/), 




e = — 



e = — 



(0.01425 
(0.01087 



iV/0.00196) GeV 4 , 
iV/0.00150) GeV\ 




(1.3) 
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where, obviously, the first and second values are due to upper and lower bounds in (1.1), 
respectively. Let us remind that these numerical values have been obtained by approximating 
the full gluon propagator by its non-perturbative term in the whole range, i. e. it has been 
already assumed that the perturbative contribution has been already subtracted. Let us 
recall that here and further on below Nf is the number of light flavors and in what follows 
we will use e = e g = eyu in (1.2-1.3) in the case of pure Yang-Mills (YM) fields, Nf = 0. 

On the other hand, many models of the QCD vacuum involve some extra classical color 
field configurations (such as randomly oriented domains of constant color magnetic fields, 
background gauge fields, averaged over spin and color, stochastic colored background fields, 
etc) and ingredients such as color-magnetic and Abelian-projected monopoles (see Refs. 
[1,5] and references therein). The relevance of center vortices to QCD by both lattice [6] 
and analytical method [7] was recently investigated as well. However, the most elaborated 
classical models are the random and interacting instanton liquid models (RILM and IILM) 
of the QCD vacuum [8]. They are based on the existence of the topologically nontrivial 
instanton-type fluctuations of gluon fields, which are solutions to the classical equations of 
motion in Euclidean space [8] (and references therein). 

In this paper we treat the chiral QCD vacuum as consisiting mainly of the two com- 
ponents, classical one given by RILM [8] and quantum one given by ZME [4] by summing 
up their contributions into the total, non-perturbative vacuum energy density. The main 
purpose of this paper is to show that this model of the non-perturbative QCD vacuum is in 
fair agreement with phenomenology. For example, it exactly reproduces the phenomenolog- 
ical value of the topological susceptibility. In sections 2, 3 and 4 using the trace anomaly 
relation [9], low-energy theorem [10,11] and Witten-Veneziano (WV) formula [12] we develop 
an analytical formalism which allows us to calculate the gluon condensate, the topological 
susceptibility and the mass of the rf meson in the chiral limit as functions of the total, non- 
perturbative vacuum energy density (the bag constant, apart from the sign, by definition). 
In section 5 we present our estimate of the non-perturbative vacuum energy density in the 
chiral limit due to instantons. Section 6 is devoted to discussion and our conclusions are 
given in section 7. The numerical results are shown in Tables 1-9. 

II. THE GLUON CONDENSATE IN THE STRONG COUPLING LIMIT 

The vacuum energy density is important in its own right as the main characteristic of the 
non-perturbative vacuum of QCD. Furthermore it assists in estimating such an important 
phenomenological parameter as the gluon condensate, introduced in the QCD sum rules 
approach to resonance physics [13]. The famous trace anomaly relation [9] in the general 
case (nonzero current quark masses m°) is 



where 6^ is the trace of the energy-momentum tensor and G a ^ u being the gluon field strength 
tensor while a s = g 2 /47r. Sandwiching (2.1) between vacuum states and on account of the 
obvious relation (0|© w |0) = 4e t , one obtains 




s 



f 



(2.1) 
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<* = ^0|^G^|0> + ^E^/(OlwlO), (2-2) 

where e t is the sum of all possible independent, non-perturbative contributions to the vacuum 
energy density (the total vacuum energy density) and (0\q^qf\0) is the quark condensate. 
From this equation in the chiral limit (m° = 0), one obtains 



(G 2 ) ee "(^G 2 ) == _(0|^G^|0> = -4e t , (2.3) 



where we need to introduce a new quantity, namely the gluon condensate in the strong 
coupling limit, i. e. not using in general the weak coupling limit solution to the /9-function 
(see section 5 below). If confinement happens then the /5-function is always in the domain 
of attraction (i. e. always negative) without IR stable fixed point [14]. Thus the non- 
perturbative gluon condensate (G 2 ), defined in (2.3), is always positive as it should be. 
Saturating e t by our values (1.2-1.3) which are relevant in the strong coupling limit, one 
obtains 

G 2 = (G 2 ) = -A{e g + N f e q ), (2.4) 

which gives the gluon condensate in the strong coupling limit as a function of Nf. The 
numerical results are shown in Table 1. 



III. THE TOPOLOGICAL SUSCEPTIBILITY 

One of the main characteristics of the QCD non-perturbative vacuum is the topological 
density operator (topological susceptibility) in gluodynamics (Nf = 0) [12] 

Xt = Jmi / #xe*±(0\T{q(x)q{0)}\0), (3.1) 

where q(x) is the topological charge density, defined as q(x) = (a s / 1 Ati)F(x)F(x) = 
(ct s /47r)F" {x)F* (x) and F^(x) = (l/2)e> J - upa Fp cr (x) is the dual gluon field strength tensor, 
while N c is the number of different colors. In the definition of the topological susceptibility 
(3.1) it is assumed that the corresponding regularization and subtraction of the perturbative 
contribution have been already done in order (3.1) to stand for the renormalized, finite and 
the non-perturbative topological susceptibility (see Refs. [10-12,15]). Precisely this quantity 
measures the fluctuation of the topological charge in the non-perturbative vacuum. 
The anomaly equation in the WV notations is 

d,J» = N f (2/N c )Q)FF. (3.2) 

As it was shown in Refs. [10,11], the topological susceptibility can be related to the 
non-perturbative gluon condensate via the low energy theorem in gluodynamics as follows 

lhru f d 4 xe^(0\T{^GG(x)^GG(0)}\0) = -£ 2 <^G 2 >. (3.3) 
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There exist two proposals to fix the numerical value of the coefficient £. The value £ = 
2/6, b = 11 was suggested long time ago by Novikov, Schifman, Vanshtein and Zakharov 
(NSVZ), who used the dominance of self-dual fields hypothesis in the YM vacuum [10]. A 
second one, £ = 4/36, was advocated very recently by Halperin and Zhitnitsky (HZ), using a 
one-loop connection between the conformal and axial anomalies in the theory with auxiliary 
heavy fermions [11] (and references therein). However, in our numerical calculations we will 
use both values for the £ parameter since the precise validity neither of the WV formula 
(4.3) below nor the NSVZ-HZ low energy relation (3.3) is not completely clear. 
The anomaly equation in the NSVZ-HZ notations is 

= N f (^)GG, (3.4) 

with Nf = 3. Thus in order to get the topological susceptibility in the WV form from 
the relation (3.3), it is necessary to make a replacement in its left hand side as follows: 
GG — > (2/N c )FF in accordance with anomaly equations (3.2) and (3.4). Then the WV 
topological susceptibility (3.1) finally becomes 

Xt = -e(^G*) = -(20 2 e YM , (3.5) 

where the second equality comes from Eqs. (2.3) and (2.4) at Nf = 0. The significance of this 
formula is that it gives the topological susceptibility as a function of the non-perturbative 
vacuum energy density for pure gluodynamics, ey hi- 
lt is easy to show that one obtains the same expression for the topological susceptibility 
(3.5) as a function of the non-perturbative YM vacuum energy density if one would use the 
weak coupling limit solution to the /3-function from the very beginning (see section 5 below). 
The numerical results due to ZME model are shown in Table 2. In conclusion, let us note 
that there exists an obvious relation between the HZ and NSVZ values for the £ parameter, 
namely iuz = (2/3)£ NS vz- 

IV. THE 17(1) PROBLEM 

The topological susceptibility (3.1) assists in the resolution of the U(l) problem^ [16] via 
the WV formula for the mass of the rj meson [12]. Within our notations it is expressed (in 
the chiral limit ) as follows: f^rni, = ANfXt, where f v > is the rf residue defined in general 

as (0\E q = u ,d,s QWM) = ijNffifP* and (0\N f ^FF\ V f ) = (iV c ^/2)/^mJ [12]. Using 
also the normalization relation f v i = \/2F®, one finally obtains 

Flm% = 2N fXt . (4.1) 

Then Eq. (3.5) implies 



1 We are going to consider here only one aspect of this problem, namely the large mass of the rf 
meson. 
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m\ = -2N f (—) e YM . 



(4.2) 



In previous expressions we omit for simplicity the superscript "0" in the pion decay constant 
as well as in m 2 ,. In the numerical evaluation of the expression (4.2), we will put, of course, 
Nf = N c = 3, while the topological susceptibility will be evaluated at Nf = as it should 
be by definition. This equation expresses the mass of the rj' meson as a function of the 
non-perturbative vacuum energy density which allows one to easily calculate it in the chiral 
limit within our formalism (see again Table 2). 

It is instructive to reproduce the WV formula (4.1) in the non-chiral case as well, namely 

m% = —f X t + A, (4.3) 

where A = 2m 2 K — m 2 . The precise validity of the WV formula (4.3) is, of course, not com- 
pletely clear, nevertheless, let us regard it (for simplicity) as exact . Using now experimental 
values of all physical quantities entering this formula, one obtains that the phenomenological 
("experimental") value of the topological susceptibility is 

^phen = a0010 58 GeV A = (180.36 MeV) 4 . (4.4) 

In the chiral limit A = since K ± and r\ particles are Nambu-Goldstone (NG) bosons. 
Omitting formally this contribution from the right hand side of Eq. (4.3) and on account of 
(1.1), one is able to derive an upper and absolute lower bounds for the mass of the r/ meson 
in the chiral limit 

854 < m\, < 913.77( MeV), (4.5) 

which should be compared with its experimental value m^f p = 957.77 MeV. One can 
conclude that the mass of rj' meson remains large even in the chiral limit. It is worth noting 
that neither the numerical value of the topological susceptibility nor the mass of the rf meson 
in the chiral limit can not exceed their phenomenological and experimental values. So the 
WV formula (4.3) in the chiral limit provides an absolute lower bound for the pion decay 
constant in this case, namely F° > 83.2 MeV. 

In order to directly apply this formalism to RILM we need the realistic estimate of the 
corresponding chiral vacuum energy density in this model. 



V. THE VACUUM ENERGY DENSITY DUE TO INSTANTONS 

The instanton-type topological fluctuations, being a classical phenomena, nevertheless 
also contributed to the vacuum energy density through a tunneling effect which was known to 
lower the energy of the ground state [8]. It can be estimated as follows. Let us consider again 
the trace anomaly relation (2.2) in the chiral limit, i. e. et = (l/4)(0|(/3(a s )/4a s )G^,G^jO). 
Using the weak coupling limit solution to the /^-function now 

a 2 2 
P(a.) = -b^ + 0(afi, b = 11 --N f , (5.1) 



6 



one obtains 

^ = -^xi(0|^G«jO). (5.2) 
The phenomenological analysis of QCD sum rules [13] for the gluon condensate implies 

G 2 = (G 2 ) = (-G 2 ) = (0|^G^J0) ~ 0.012 GeV 4 , (5.3) 

71 71 

which can be changed within a factor of two [13]. From the phenomenological estimate (5.3), 
one easily can calculate (l/8)(0|(o; s /7r)G'^G'^|0) ~ 0.0015 GeV 4 ~ 1.0 jmr 4 . Having in 
mind this and assuming that the gluon condensate in the weak coupling limit is deter- 
mined by the instanton-type fluctuations only, Shuryak [8] (see also references therein) has 
concluded in that the "average separation" between instantons was R ~ 1.0 fm, so the cor- 
responding density of the instanton-type fluctuations should be n ~ 1.0 fm~ 4 . Let us note 
that the second parameter of the instanton liquid model of the QCD vacuum, the instanton 
size po — 1/3, was chosen to reproduce standard (also (as gluon condensate) phenomenolog- 
ically estimated from QCD sum rules [13]) value of the quark condensate. This contribution 
to the vacuum energy density via the trace anomaly relation (2.1-2.2) vanishes in the chiral 
limit. However, due to all reasonable estimates of light quark masses, numerically its con- 
tribution is at 20% and thus comparable with the systematic error in the determination of 
the gluon condensate itself [13,17]. 

Saturating the total vacuum energy in the weak coupling limit by instanton component ej 
and using the above mentioned estimate, from (5.2) for dilute ensemble, one finally obtains 

e, = ~n = ~ x 1.0 fm~ 4 = -(0.00417 - iV/0.00025) GeV 4 . (5.4) 

Thus instanton contribution to the vacuum energy density was not calculated independently 
but was postulated via the trace anomaly relation using the phenomenological value of the 
gluon condensate (5.3) as well as weak coupling limit solution to the /5-function (5.1). It is 
well known that density of instanton-type fluctuations is suppressed in the chiral limit and 
is again restored bacause of dynamical breakdown of chiral symmetry [8] (and references 
therein). In any case it can not be large in the chiral limit, so the functional dependence of 
the vacuum energy density on the instanton density, established in Eq. (5.4) due to dilute 
gas approximation, seems to be justified in this case. The only problem is the numerical 
value of the instanton density itself, which can be taken either from phenomenology or from 
lattice simulations. 

In Ref. [10] it was argued that the gluon condensate in the chiral limit is approximately 
two times less than the above mentioned phenomenological (empirical) value (5.3), i. e. 
(G 2 ) ch ~ 0.5(G 2 ) phen . This means that in this case instanton density n ~ 0.5 fm~ 4 and the 
vacuum energy density due to instantons approximately two times less than (5.4). How- 
ever, it has been already pointed out [18] that QCD sum rules substantially underestimate 
the value of the gluon condensate. The most recent phenomenological calculation of the 
gluon condensate is given by Narison in Ref. [19], where a brief review of many previous 
calculations is also presented. His analysis leads to the update average value as 

(0|— G a G a |0) = (0.0226 ± 0.0029) GeV 4 . (5.5) 
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This means that instanton density is approximately two time bigger than it was estimated 
by Shuryak for instanton liquid model [8], but in the chiral limit we are again left with (5.4). 

In lattice QCD situation with instaton density and their sizes is also ambigious. In 
quenched (Nf = 0) lattice QCD by using the so-called "cooling" method the role of the 
instanton-type fluctuations in the QCD vacuum was investigated [20]. In particular, it 
was found that the instanton density should be n = (1 + 5) fm~ A , where 5 ~ 0.3 — 0.6 
depending on cooling steps. Moreover, by studying the topological content of the vacuum 
of 577(2) pure gauge theory using a method of RG mapping [21], it is concluded that the 
average radius of an instanton is about 0.2 fm, at a density of about 2 fm~ A . However, 
in Ref. [22] the topological content of the SU(3) vacuum was studied using the same 
method as for SU(2) gauge theory earlier and was obtained a fair agreement with Shuryak's 
phenomenologically estimated numbers for the instanton liquid model. At the same time, 
in Refs. [23,24] considerably larger values were reported. Thus at this stage it is rather 
difficult to choose some well-justified numerical value of the instanton-type contribution to 
the non-perturbative vacuum energy density. In any case, in what follows we will consider 
(5.4) as a realistic upper bound for the instanton contribution to the vacuum energy density 
in the chiral limit. If the instanton number density is about n ~ 2 fm~ A , then in the chiral 
limit we again are left with (5.4), but if it is about n ~ 1 fm~ A , we will be left with half of 
(5.4). Then the instanton contributions to the topological susceptibility and the mass of the 
rj' meson in the chiral limit are to be calculated via (3.4) and (4.2), respectively on account of 
the substitution tyu — ► £i(Nf = 0), where e/ is given in (5.4) with the two different values 
for instanton number densities in the chiral limit, n = 0.5 fm~ A , 1.0 fm~ A . The numerical 
results are shown in Table 3. In conclusion, we note that for densities n > 2 fvrC A (which 
means n > 1 fvrT A in the chiral limit) the applicability of the dilute gas approximation 
becomes, apparently, doubtful. 



VI. DISCUSSION 

A. The gluon condensate 

It becomes almost obvious that we must distinguish the two types of gluon condensates, 
both of which are the non-perturbative quantities. The first one is determined by (2.3) and 
is the one which is relevant in the strong coupling limit. In this case the total vacuum energy 
is mainly saturated by the ZME component as it is precisely shown in (2.4). In the weak 
coupling limit, saturating e t by ej, from (5.2-5.4) one obtains 

a s 2 32 32 

( — G ) = — -e t ~ — -e/ = 8n, (6.1) 

71 

i. e. gluon condensate in the weak coupling limit does not explicitly depend on Nf. As 
was mentioned above, precisely this gluon condensate was introduced long ago [13]. This 
unphysical situation takes place because in instanton calculus [8] there is no other way 
to calculate the vacuum energy density than the trace anomaly relation (2.1-2.2) which 
becomes finally (6.1) as it was described above. In this case it is preferable to have the Nf 
dependent vacuum energy density than the gluon condensate since the former is the main 
characteristic of the non-perturbative vacuum. Contrast to this, we have calculated the 
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vacuum energy density completely independently from the trace anomaly relation. We use 
it only to calculate the gluon condensate in the strong coupling limit. That is why in our 
case both quantities are Nf dependent functions. 

Let us make a few things perfectly clear. It makes sense to underline once more that 
the vacuum energy density is not determined by the trace anomaly relation (2.2). The real 
situation is completely opposite. As it was already mentioned in the Introduction, the QCD 
vacuum beeing a very complicated medium, may contain many different components and 
ingredients which contribute to the vacuum energy density. These contributions are com- 
pletely independent from the gluon condensate, of course. For example, in the chiral limit 
the explicit quark contribution to the vacuum emergy density via the trace anomaly rela- 
tion (2.2) vanishes . However, because of the dynamical chiral symmetry breaking (DCSB), 
there is nonvanishing (even in the chiral limit) explicit contribution to the quantum part of 
the vacuum energy density which comes from the vacuum quark loops as it was described 
in our work [4] in detail. Thus, the total vacuum energy density, defined as the trace of 
the energy-momentum tensor, becomes the sum of all independent contributions. This sum 
precisely determines the realistic value of the gluon condensate in the chiral limit via the 
trace anomaly relation and not vice versa. The gluon condensate may exist or not, but the 
vacuum energy density always exists (at least its quantum part) due to non-perturbative 
solutions to quark and gluon SD equations by substituting them into the effective potential 
(see Ref. [4] again). In other words, the vacuum energy density is much more fundamental 
quantity than the gluon condensate. 

Our bounds for full QCD (Nf = 3) gluon condensate in the strong coupling limit 



are comparable with recent phenomenological determination of the standard gluon conden- 
sate by Narison (5.5). The parameterization (the left hand side) of the two types of the 
gluon condensate may be, of course, the same but their numerical values (the right hand 
sides) are not to be the same. This difference is not only due to different physical observables 
as was noticed in Ref. [8]. Though both quantities are the non-perturbative phenomena, 
nevertheless this difference reflects different underlying physics. Our gluon condensate (2.3) 
is the strong coupling limit result and reflects the nontrivial topology of the true QCD vac- 
uum where quantum excitations of gluon fields play an important role. As was shown in our 
preceding papers [4,25] precisely these type of gluon field configurations are mainly respon- 
sible for quark confinement and DCSB. At the same time, the standard gluon condensate 
(5.3) is the weak coupling phenomenon due to classical instanton-type fluctuations in the 
true QCD vacuum which by themselves do not confine quarks [21,26-28]. 

Concluding let us note that in the lattice simulations there already exist calculations 
of the gluon condensate which are one order of magnitude bigger than the standard value, 
namely G 2 ~ 0.1046 GeV A for SU(3) in Ref. [29] and G 2 ~ 0.1556 GeV 4 for SU(2) in Ref. 
[30] ( see also review [8]). In phenomenology also there exist large values, namely 



0.025 < -(0| 



f3(a s ) 



G* GIM < 0.033 (GeV 4 ), 



(6.2) 



4a s 



0.04 < (OpG; G* |0) < 0.105( GeV 4 ), 



(6.3) 



which were recently derived from the families of J/ty and T mesons in Ref. [31]. 
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B. Topology of chiral QCD vacuum 



Our numerical results for the quantum part of the non-perturbative vacuum energy 
density and for the topological susceptibility with the mass of the 7/meson in the chiral limit 
are presented in Eqs. (1.2-1.3) and in Table 2, respectively. In general our values for the 
vacuum energy density are an order of magnitude bigger than RILM can provide at all in 
various modifications (see Table 3). That is why the quantum part of the non-perturbative 
vacuum energy density saturates the phenomenological value of the topological susceptibility 
and the mass of the rj' meson in the chiral limit much better than the classical part given 
by instantons (compare Tables 2 and 3). Especially this is obvious for the HZ value of the 
£ parameter, introduced in the low-energy theorem, Eq. (3.3). The instanton contribution 
substantially underestimates the phenomenological value of the topological susceptibility 
and therefore can not account for the large mass of the rj' meson in the chiral limit alone 
(see Table 3). 

However, the total vacuum energy density, e t , is, in principle, the sum of all possible 
independent, the non-perturbative contributions. Thus, at least it is the sum of the two 
well-established contributions, quantum e = tzME and classical ej, i. e. et — €zme + £1 + 
where the dots denote other possible independent contributions. In this case an excellent 
agreement with phenomenology is achieved indeed (see Tables 4 and 5). The numerical 
values of the bag constant B, defined as the difference between the perturbative and non- 
perturbative vacua are given now by the relation B = —e t and can be explicitly evaluated 
using Eqs. (1.2-1.3) for €zme and Eq. (5.4) for ej on account of the above mentioned two 
different instanton number densities (see Tables 6-9). For the readers convenience the bag 
constant (and consequently the total, non-perturbative vacuum energy density) is given in 
often used different physical units. 

VII. CONCLUSIONS 

In summary, using the trace anomaly relation, NSVZ and HZ low-energy theorem and 
Witten-Veneziano formula, we have developed an analytical formalism which allows one to 
calculate the gluon condensate, the topological susceptibility and the mass of the rj' meson 
in the chiral limit as functions of the non-perturbative vacuum energy density. It was 
immediately used for numerical investigation of the chiral QCD non-perturbative vacuum 
topology within the recently proposed ZME quantum model. We have explicitly shown 
that precisely our values for the non-perturbative vacuum energy density (1.2-1.3) are of 
the necessary order of magnitude in order to saturate the large mass of the rj' meson in 
the chiral limit. We have obtained good approximation to the phenomenological value of 
the topological susceptibility as well. The HZ value of the £ parameter, introduced in the 
low-energy theorem (3.3), especially nicely saturates them (for all result mentioned above 
see Table 2). At the same time, it is clear that instanton-induced contribution should be 
added to our values in order to achive an excellent agreement with phenomenology. Indeed, 
from Table 5 it follows 

Xt = (180.3 MeV f (7.1) 
and consequently (as it should be) 
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m\, = 854 MeV, 



(7.2) 



which are in fair agreement with (4.4) and lower bound in (4.5), respectively. Let us re- 
mind that these numbers have been obtained when the pion decay constant was precisely 
approximated by its experimental value. The above displayed excellent agreement with 
phenomenological values of the corresponding quantities is achieved by summing up our 
contribution and instanton-induced contribution into the total, non-perturbative vacuum 
energy density, i. e. the summation was done purely phenomenologically by simply sum- 
ming up the two well-established contributions (classical instanton's and quantum ZME's). 
How to take into account 't Hooft's instanton-induced interaction [32] at the fundamental 
quark level within our approach is not completely clear for us, though see paper [33] (and 
references therein). 

Let us also make a few things perfectly clear. At low energies QCD is governed by 
SU^Nf) x SUn(Nf) chiral symmetry and its dynamical breakdown in the vacuum to the 
corresponding vectorial subgroup [34]. The chiral limit is not physical one but nevertheless 
remains very important theoretical limit since to understand the chiral limit physics means 
to correctly understand the dynamical structure of low-energy QCD as well as the topolog- 
ical properties of its ground state. So a realistic calculation of various physical quantities 
as well as chiral properties of its vacuum becomes important. In particular, any model 
of the QCD vacuum should pass the chiral limit test in order to be justified for further 
extrapolation to the realistic (nonchiral) case. In our previous publication [4] ZME quan- 
tum model was formulated. Here we have explicitly shown its important and novel feature, 
namely it itself passes the chiral limit test justified thereby for use in the nonchiral case 
as well. Complemented by instanton-induced contribution it is in a fair agreement with 
phenomenology. 

In conclusion a few remarks are in order. It is well-known that instanton-type fluctuations 
require topological charge to be integer (±1) and the vacuum angle, 9 - nonzero, which 
violates P and CP invariance of strong interactions [10,16]. The non-perturbative q~ A - 
type quantum excitations do not require the introduction of the vacuum angle, 9, at all. 
It is quite possible that topological charge in this case is not restricted to integer values. 
It has been explicitly shown that fractional (non-integer) topological charge configurations 
are required to resolve the U(l) problem [16,35]. However, the 9 dependence of the QCD 
non-perturbative vacuum energy remains an important problem. For recent developments 
of this problem in the large N c limit of four-dimensional gauge theories see papers [36]. In 
particular, in Refs. [37,38] it has been discussed that the picture of its dependence in QCD 
for finite N c might be more complicated than that predicted by the large N c values. 

And finally, let us emphasize once more indisputable simplicity of our analytical calcu- 
lation of the topological susceptibility (7.1) in comparison with indisputable complexity of 
its calculation by lattice method [21-24,30,39-41]. 
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Table 1: The gluon condensate in the strong coupling limit 



F° (MeV) G 2 (GeV 4 ) 


N f = 


N f = l 


N f = 2 


N f = 3 


87.2 


0.043 


0.037 


0.031 


0.025 


93.3 


0.057 


0.049 


0.041 


0.033 



Table 2: ZME model values for Xt an d rn v i in MeV units 



F° 

7T 


Xt" (NSVZ) 


(HZ) m°, 


(NSVZ) 


(HZ) 


87.2 


195.4 


159 


1065 


710 


93.3 


208 


170.24 


1139 


759.7 



1 



Table 3: Chiral topology due to instantons 



n (fm 4 ) 


0.5 


1.0 


ei (GeV 4 ) 


-0.002085 


-0.00417 


G 2 (GeV 4 ) 


0.006 


0.012 


1 1 A 

Xt V4 (MeV) 






(NSVZ) 


128 


153 


(HZ) 


104 


125 


ra°, (MeV) 






F° = 87.2 MeV 






(NSVZ) 


466.5 


659.55 


(HZ) 


311 


439.7 


m°, (MeV) 






F° = 93.3 MeF 






(NSVZ) 


435.9 


616.5 


(HZ) 


290.6 


411 



2 



Table 4: Chiral QCD topology (F° = 87.2 MeV) 





£zme + e/(0.5 /m 4 ) 


ezME + e/(1.0 fm 4 ) 


e t (GeV 4 ) 


-0.013 


-0.015 


Xt V4 (MeV) 






(NSVZ) 


203.5 


211 


(HZ) 


166.2 


172.3 


ra°, (MeV) 






(NSVZ) 


1164.6 


1251 


(HZ) 


776.4 


834 



3 



Table 5: Chiral QCD topology (F° = 93.3 MeV) 



£zme + ej(0.5 /m 4 ) £zme + e/(1.0 /m 4 ) 



e* (GeV 4 ) 


-0.016 


-0.018 


Xt V4 (MeV) 






(NSVZ) 


214 


220.8 


(HZ) 


175 


180.3 


raJJ, (MeV) 






(NSVZ) 


1207.5 


1281 


(HZ) 


805 


854 



Table 6: The bag constant (F° = 87.2 MeV, n = 0.5 fm~ 4 ) 

B = —e t N f = N f = 1 Nf = 2N f = 3 

GeV A 0.013 0.01133 0.097 0.008 

MeV 4 (337.7) 4 (326.25) 4 (313.8) 4 (300) 4 

GeV/fm 3 1.7 1.47 1.26 1.04 



4 



Table 7: The bag constant (F° = 93.3 MeV, n = 0.5 fm^) 



B = —e t N f = N f = l N f = 2 N f = 3 

GeV 4 0.016 0.01426 0.01216 0.01 

MeV 4 (355.65) 4 (345.56) 4 (332) 4 (316) 4 

GeV/fm 3 2.0 1.8 1.58 1.30 



Table 8: The bag constant (F° = 87.2 MeV, n = 1 fm~ 4 ) 

B = —e t N f = N f = l N f = 2 N f = 3 

GeV 4 0.015 0.0133 0.01154 0.0098 

MeV 4 (350) 4 (339.6) 4 (327.75) 4 (314.6) 4 

GeV/fm 3 1.95 1.73 1.5 1.27 
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Table 9: The bag constant (F° = 93.3 MeV, n = 1 fm~ A ) 



B = —e t N f = N f = l N f = 2 N f = 3 

GeV 4 0.01842 0.0162 0.014 0.0118 

MeV 4 (368.4) 4 (356.76) 4 (344) 4 (330) 4 

GeV/fm 3 2.4 2.1 1.82 1.53 
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